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Preface: How many permutations are there of length googol+30 avoiding an increasing 
subsequence of length googol? 

This number is way too big for our physical universe, but the number of permutations of length 
googol+30 that contain at least one increasing subsequence of length googol is a certain integer that 
may be viewed in http://www.math.rutgers.edu/~zeilberg/mamarim/mamarimPDF/gessel64.pdf. 
Hence the number of permutations of length googol+30 avoiding an increasing subsequence of length 
googol is ( googol + 30)! minus the above small number. 

Counting the “Bad Guys” 

Recall that thanks to Robinson-Schensted ([Rob] [Sc]), the number of permutations of length n that 
do not contain an increasing subsequence of length d is given by 

Gd(n) := X /a , 

Ahn 

Jfcrows(X') <d 


where A denotes a typical Young diagram , and f\ is the number of Standard Young tableaux whose 
shape is A. 

Hence the number of permutations of length n that do contain an increasing subsequence of length 
d is 

B d (n) := X A 2 • 

Ahn 

#rows(\)>d 

Since the total number of permutations of length n is n\ ([B]), if we know how to find B d (n ), we 
would know immediately G d {n) = n\ — B d (n ), at least if we leave n! alone as a factorial, rather 
than spell it out. 

Recall that the Hook Length formula (see [Wiki]) tells you that if A is a Young diagram then 


Ucex Me) 

where the product is over all the n cells of the Young diagram, and the hook-lenght , h(c), of a cell 
c = (i,j), is (Ai — i) + (A' — j) + 1, where A' is the conjugate diagram, where the rows become 
columns and vice-versa. 
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Let r be a fixed integer, then for symbolic d, valid for d > r — 1, any Young diagram with at least 
d rows, and with d + r cells, can be written, for some Young diagram p = (/q,..., fx r ), with < r 
cells, (where we add zeros to the end if the number of parts of /j is less than ?’) as 

^ = (1 + Hi, . . • , 1 + Hr, \ d ~ r+r ') , 

where r' = r — |/x|. For such a shape A, with at least d rows, 


II h ( c ) = ( II 

c£A \ c£/.i 

Hence /a, that is (d+r)! divided by the above, is a certain specific number times a certain polynomial 
in d. Since, for a specific, numeric, r, there are only finitely many Young diagrams with at most r 
cells, the computer can find all of them, compute the polynomial corresponding to each of them, 
square it, and add-up all these terms, getting an explicit polynomial expression, in the variable d, 
for B d (d + r), the number of permutations of length d + r that contain an increasing subsequence 
of length d. As we said above, from this we can find G d {d + r) = (d + r)! — B d (d + r ), valid for 
symbolic d > r — 1. 


^ • ((d + r' + Hi)(d + r' - 1 + // 2 ) ■ ■ ■ (d + r' - r + 1 + n r )) ■ (d - r + r')\ 


Bd(d + r) for r from 0 to 


3U 


B d (d) = 1 , 

B d (d + 1) = d 2 + 1 
1 


B d (d + 2) — — d 4 + d 3 + — d 2 + d + 3 , 


2 2 
5 .a 2 


B 


B d (d + 3) = - d 6 + d 5 + - d 4 + - d 3 + — d 2 + — d + 11 , 

y ’ 6 3 3 6 3 

1 


1 

2 

B d (d+5 ) = — d 10 +i d 9 +- d 8 +- d 7 +— d 6 +- d 5 d 
v ; 120 6 24 3 120 30 24 

For B d (d + r) for r from 6 up to 30, see 

http://www.math.rutgers.edu/~zeilberg/tokhniot/oGessel64a . 
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Sequences 


The sequence Gs(n) is the greatest celeb in the kingdom of combinatorial sequences [the subject 
of an entire book([St]) by Ira Gessel’s illustrious academic father, Richard Stanley], the super- 
famous A000108 in Neil Sloane’s legendary database ([SI]). G^{n), while not in the same league 
as the Catalan sequence, is still moderately famous, A005802. G${ri) is A047889, Gq(u) is 
A047890, G 7 (n) is A052399 , G 8 (n) is A072131, G 9 (n) is A072132, G w {n) is A072133, Gn(n) 
is A072167, but G d (n) for d > 12 are absent (for a good reason, one must stop somewhere!). Also 
the flattened version of the double-sequence, {G d {n)}, for 1 < d < n < 45 is A047887. Using the 
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polynomials Bd(d + r), we computed the first 2d + 1 terms of Gd(n ) for d < 30. See 

http://www.math.rutgers.edu/~zeilberg/tokhniot/oGessel64b . 

But this method can only go up to 2d + 1 terms of the sequence G^(n), and of course, the first 
d — 1 terms are trivial, namely n\. Can we find the first 100 terms (or whatever) for the sequences 
Gd(n) for d up to 20, and beyond, efficiently? 

Encore: Efficient Computer-Algebra Implementation of Ira Gessel’s AMAZING De¬ 
terminant Formula 


Recall Ira Gessel’s [G] famous expression for the generating function of Gd(n)/n\ 2 , canonized in the 
bible ([W], p. 996, Eq. (5)). Here it is: 

H x2n = det ( J |i-il ( 2x ))i,j =i.rf i 

n> 0 


in which I v {t) is (the modified Bessel function) 


«)-£r 


j-r'o iKd + »)'■ 


Can we use this to compute the first 100 terms of, say, G 2 o{n)l 


While computing numerical determinants is very fast, computing symbolic ones is a different story. 
First, do not get scared by the “infinite” power series. If we are only interested in the first N terms 
of Gd(n), then it is safe to truncate the series up to t 2N , and take the determinant of a dx d matrix 
with polynomial entries. If you use the vanilla determinant in a computer-algebra system such as 
Maple, it would be very inefficient, since the degree of the determinant is much larger than 2 N. 
But a little cleverness can make things more efficient. The Maple package Gessel64, available free 
of charge from 


http://www.math.rutgers.edu/~zeilberg/tokhniot/Gessel64 , 

accompanying this article, has a procedure SeqIra(k,N) that computes the first N terms of Gk{n), 
using a division-free algorithm (see [Rot]) over an appropriate ring to compute the determinant in 
Gessel’s famous formula. 


SeqIra:=proc(k,N) local ira,t,i,j, R: 
R := table(): 

R[‘0‘] := 0: 

R[‘1‘] := 1: 

R[‘+‘] ;= 
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R [ ! -‘] := : 

R[ £ * c ] := proc(p, q): return add(coeff(p*q, t, i)*t**i, i=0..2*N): end: 

R[ £ = < ] := proc(p, q): return evalb(p = q): end: 

ira:=expand(LinearAlgebra[Generic][Determinant][R](Matrix([seq([seq(Iv(abs(i-j),t,2*N), 

j=l..k-1)], 
i=l..k-l)]))): 


[seq(coeff(ira,t,2*i)*i!**2,i=l..N)] : 
end: 


In the above code, procedure Iv(v,t,N) computes the truncated modified Bessel function that 
shows up in Gessel’s determinant, and it is short enough to reproduce here: 

Iv:=proc(v,t,N) local j: add(t**(2*j+v)/j!/(j+v)!,j=0..trunc((N-v)/2)+l): end: 

Using this procedure, the first-named author computed (in 4507 seconds) the first 100 terms of 
each of the sequences Gd(n ) for 3 < d < 20, and could have gone much further. 

See http://www.math.rutgers.edu/~zeilberg/tokhniot/oGessel64c 

HAPPY 64th BIRTHDAY, IRA! 
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